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Abstract 

We use NRQCD to predict the cross section for r] c production at the Tevatron. 
The non-perturbative matrix elements required are obtained, using heavy-quark 
symmetry, from the matrix elements determined from the CDF J /ip data. Our 
numbers are, therefore, predictions of NRQCD and provide a very good test of 
it. Even after taking into account the small branching ratio of the n c into two 
photons, a substantial rate for the production of this resonance at the Tevatron 
is expected. 
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It is by now well known from compelling theoretical QTJ| and phenomenological 
reasons that the Colour Singlet model (CSM) used to describe the production and 
decay of quarkonium is incomplete. A new factorisation approach within QCD was 
provided by Bodwin, Braaten and Lepage ||. This method separates the perturba- 
tively calculable short- distance (< 1/Mq, Mq is the mass of the heavy quark) effects 
from the long-distance effects, which are described by non-perturbative matrix ele- 
ments in an effective field theory called non-relativistic QCD (NRQCD). Unlike CSM 
where factorisation is explicitly violated, factorisation can be shown to be valid in 
NRQCD, because in this approach the bound state problem is treated in a system- 
atic expansion in powers of v, the relative velocity between the quarks in the bound 
state. Consequently, in NRQCD the production and decay of quarkonium takes place 
via intermediate states in which the QQ pair has quantum numbers different from 
those of the physical quarkonium. The QQ pair in the intermediate state denoted as 
QQ[ 2S+1 L^' ] could be in the colour singlet or octet state and could form the physical 
quarkonium state by making chromo-electric or -magnetic transitions. Thus the cross 
section for production of a quarkonium state H can be written as 
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The coefficients F n correspond to the production of QQ in the angular momentum 
and colour state denoted by n and is calculated using perturbative QCD. The non- 
perturbative part, (O^) of mass dimension d n in NRQCD has a well-defined operator 
definition and is universal. These matrix elements can be extracted from any one 
process and then be used to predict other processes where the same matrix elements 
appear. Though the summation involves an infinite number of terms, the relative 
magnitude of the various terms is predicted by NRQCD as these matrix elements 
scale as powers of v. However, this does not necessarily imply that effects from higher 
orders in v will always be small in physical processes, because any observable, like the 
decay width or the cross section is given by a double expansion in the strong coupling 
constant a s (Mo) and the relative velocity v. Indeed, the higher order colour-octet 
components have important effects in certain processes: for example, in \ c decays 
into light hadrons, where these octet components are necessary to obtain an infra- 
red safe expression for the decay-width ||, and in large-p^ J/i> and ip' production 
from the CDF experiment at Tevatron where discrepancies larger than one order 
of magnitude vis a vis CSM predictions could be explained by including the octet 
components [|]. 0. 

The non-perturbative matrix elements in NRQCD, are not calculable and have to 
be obtained by fitting to available data. The matrix elements of the colour-singlet 
operators are related to the radial wave-functions and can be obtained from decay 
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widths or lattice calculations, for example, but the colour-octet matrix elements are 
not known and have only been obtained by fitting NRQCD predictions to the CDF 
data H 0| . Given that the normalisation of the cross sections are not predictions 
of the theory, it is important to look for other independent quantitative tests of 
NRQCD. One crucial test of NRQCD was thought to be the production of large-pr 
inelastic J /if) production in ep collisions at HERA. The inelastic cross section has 
been measured by both HI and ZEUS ||. The inelasticity of the events is ensured 
by selecting events with z < 1, where z = p p ■ pj/^/p p ■ p 7 . In addition, a px cut 
of 1 GeV is used to select the events. Using the values of the non-perturbative 
inputs obtained from the Tevatron, the theoretical predictions have been compared 
H to these inelastic data. However, there seems to be a problem in predicting the z 
distribution because while the colour-singlet cross section dominates in most of the 
\ow-z region, the colour-octet contribution grows anomalously in the large- z (0.8 < 
z < 0.9) region and this rise is not seen in the data. However, it is premature to 
conclude that this is a failure of NRQCD, because at z values of 0.8 other effects 
like soft-gluons [10], fc^-smearing |l ij could become important, or even worse, the 



NRQCD factorisation may break down . The safer conclusion is that the inelastic 



J/if> production process at HERA is not a clean test of NRQCD. 

In fact, NRQCD is best tested at Tevatron itself in other large-pr observables. 



The polarisation of the produced J/if> is one such test |T3|- The second test is to 
look for the production of other charmonium states at large-pr at the Tevatron. This 
is the philosophy of the present paper and we study \aige-p T r\ c production at the 
Tevatron, with this end in mind. The remarkable thing is that the non-perturbative 
parameters appearing in the n c production cross section can be determined from the 
matrix elements determined from J /if} production at the Tevatron: this happens 
because of the heavy-quark symmetry of the NRQCD Lagrangian. This has been 
exploited earlier in the context of h c production at Tevatron |14]]. We show, in the 
following, that from NRQCD one expects a large cross section for n c production at 
the Tevatron. The canonical search channel for the n c is its two photon decay mode 
and the number of events is substantial even after taking the small branching ratio 
of the rj c — > 77 into account. 

By writing a Fock space expansion of the physical r] c , which is a 'So (J PC = 0-+) 
state, we have 

\ Vc ) = 0(1) QW4 11 ]) + 0{v 2 ) QWP?] a) + 0(v 4 ) QQ[ 3 S?] <?) + •••■ (2) 



The colour-singlet 1 S , state contributes at 0(1) but the colour-octet x Pi and 3 Si 
channels effectively contribute at the same order because the P-state production is 
itself down by factor of 0(v 2 ). The colour-octet states become a physical r) c by the 
X P| 8 ' state emitting a gluon in an El transition, while the 3 Sj 8 ' emitting a gluon in 
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an Ml transition. To 0(a^v 7 ) the r] c production cross section is 



(0| 0?f Px] |0) + ^1 (0| Ogf |0> • (3) 

The coefficients F's are the cross sections for the production of cc pair in the respective 
angular momentum and colour states and is given by 

da { P p->cc[ 2S+l L [1 Ax) = 



dpi 

Y,J d yJ **i *i G */v(*i) G b / P {x 2 ) — g|- - ^(a6 - cc^L?' 81 ] d). (4) 
The contributing subprocess cross sections are 

- QQ[ 2S+1 Lj] g, 
gq(q) - QQ[ 2S+1 Lj] q (q), (5) 

gg - qq[ 25+ %u 

where the is in the 1 S , q 1 ', 1 p| 8 ' and 3 Sj 8 ^ states. The summation is over the partons, 
a and 6. G a /p and Gb/p are the distributions of partons a and 6 in the hadrons and 
x\ and X2 are the momentum they carry, respectively. x 2 is given in terms of x\ as 

xi x T - 2r 
^2 = , (o) 

2x x -x T e y 

where r = M 2 /s, with M the mass of the resonance, yfs the centre-of-mass energy, 



y the rapidity at which the resonance is produced and Xt = y x^ + At = 2Mt/ 
with xt = 2pr I \fs. 

Among the coefficients the matrix elements for the subprocesses corresponding to 



Fi^So] are available in fl5|] and those corresponding to F 8 [ 3 Si] are calculated by Cho 
and Leibovich 0. We have calculated the remaining coefficient Fg[ 1 Pi]. The packages 
FORM and MATHEMATICA are used to calculate the matrix element squares of the 
subprocesses. 

One important component at large-pr is the gluon fragmentation to charmonium, 
as has been shown in the case of large-pr J/ip production Q. It has been noted 
in Ref. || that even in the fixed-order calculation of the J/ip cross section, the dia- 
grams that contribute at large-pr are fragment at ion- like diagrams in which a single 
gluon attaches itself to a cc pair. In our case, we find that while the Fgp^i] has 



3 



a fragmentation-like contribution at large-pr there is no such contribution for the 
F 8 [ 1 P 1 ]. Following Ref. ||, we have also corrected the cross section by the ratio of 
Q 2 evolved to the unevolved fragmentation function to get the correct shape of the 
Pt distribution at large pr- 

Heavy quark spin-symmetry is made use of in obtaining (C^ c )'s from the exper- 
imentally available (O^'^)'s. Using this symmetry we get the following relations 
among (O^)'s: 

(o\o?[ l s ]\o) = (oio^p&iio) (i + o(v*)), 

(OlOgfPJIO) = (0|Of^P ]|0> (l + 0(v 2 )), 

(0|C»r[ 3 5i]|0) = (0|O 8 J/ ^[%]|0) (l + 0(v 2 )). (7) 

Hence, for the sing let matrix elements we have (0| O^pSi] |0) = 1.2 GeV 3 [| and 
for the octet matrix elements extracted from the CDF J/ ip data we have Ai + A 2 = 

(O|0 V [ : Pol|O) + mS!^m = (2.2 ± 0.5) x 10- 2 GeV 3 §. The CDF J/rj, do not 
allow for a separate determination of the values of A\ and A 2 because the shapes of 
these two contributions to the J / ip pt distribution are almost identical. We use the 
fact that the sum A x + A 2 is constrained by the CDF data, and for our numerical 
predictions we assume that either A\ or A 2 saturates the sum. In one case, therefore, 
we have the maximum possible contribution from the 3 Sj 8 ' channel and none from 
the x p| 8 ' channel, and in the other case the 3 S'{ 8 ' channel makes no contribution while 
the x p| 8 ' contributes its maximum. 

In Fig. 1, we have plotted the differential cross section da/dpx against pr- We are 
considering the case where the r\ c will be searched for in its decay into two photons 
and we have consequently folded in the differential cross section with the i] c —>■ 77 
branching ratio (B = 3 x 10~ 4 ). In our computations, we have used MRSD-' parton 
densities fl6|] and the parton densities are evolved to a scale Q = Mt- As a conse- 



quence of the fact that we are saturating the sum with either A\ or A 2 , our curves for 
the two octet channels represent the maximum possible cross section that can result in 
each (with the other being zero). Because of the fact that fragmentation dominates 
at high pt contributions from and 1 P| 8 ', which do not have a fragmentation 
contribution, fall steeply with increasing p?, whereas the 3 Sj^ contribution makes a 
sizable contribution even at the largest px considered. The shape of the measured pt 
distribution is clearly a very good measure of the individual matrix elements A\ and 
A 2 . We reiterate that the CDF J ftp data does not offer such a discrimination. 

To get an idea of the measurability of the rj c cross section at the Tevatron, we have 
integrated our p T distribution for values of px above 5 GeV. Assuming an integrated 
luminosity of 110 pb _1 we find that the number of events in the 77 channel from the 
singlet is about 400 and that from 3 S'{ 8 ^ is about 7300 when x p| 8 ' contribution 
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is absent, while the number of events from 1 P 1 ' 8 ' is about 25 when 3 5'{ 8 ' contribution 
is absent. Thus the number of r] c — > 77 lies between 425 and 7700. Thus, even the 
minimum of 425 events provides a very good prospect for the observation of r] c at the 
Tevatron. We also point out that the integrated cross sections are also good measures 
of the matrix elements A\ and A2. 

We have studied the effect on the cross section of the variation of the parton 
densities, the scale and the non-perturbative matrix elements. By using GRV HO 



densities | 17fl , instead of MRSD-', we find that the cross section increases by about 
25%. If we use the scale choice Q = M T /2 instead of Q = M T then the cross section 
increases by about 45-50%. In addition to these uncertainties, we can expect a 30% 
variation in the values of the non-perturbative matrix elements we have used, because 
of the fact that the heavy-quark symmetry is approximate. 

In our calculations, we have considered only the direct production of r] c . There is 
a contribution to the r] c signal coming from the decays of J/i/j's. From the measured 
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J/ip cross section at CDF we find that the contribution to the signal coming from 
J/tp decays will be roughly about 100 events which is only a small fraction of the 
minimum number of events that we expect from the direct production process. A 
more complete analysis including the J/ip decay contribution to the r\ c signal will be 
presented in a future publication. 

We conclude with the following observations: the heavy-quark symmetry of the 
NRQCD Lagrangian allows us to make predictions for the r\ c cross section at the Teva- 
tron. The integrated rate and the shape of the p T distributions are very sensitive to 
the non-perturbative matrix elements A\ and A 2 , whose sum has been determined 
from the J/ifj production data from CDF, but are not known individually. This sen- 
sitivity can clearly be used to determine these matrix elements individually. Our 
estimates for the integrated rate vary between 425 and 7700 r\ c events (after including 
the branching ratio of the rj c into two photons). Even at its minimum of 425 events, 
the prospect of observability of the r/ c signal at the Tevatron seems extremely good. 
We emphasise that this is a testable prediction of NRQCD, and that such a predic- 
tion is not possible in alternate approaches to quarkonium formation like the colour 
evaporation model. 
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